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Abstract 

Threshold neutral pion photoproduction on the deuteron is studied in the frame- 
work of baryon chiral perturbation theory beyond next-to-leading order in the chiral 
expansion. To fourth order in small momenta, the amplitude is finite and a sum 
of two- and three-body interactions with no undetermined parameters. With ac- 
curate theoretical and experimental input from the single nucleon sector for the 
proton amplitude, we investigate the sensitivity of the threshold cross section to 
the elementary 'yn — > Tr^n amplitude. A precise measurement of the threshold cross 
section for 'yd -^ n^d is called for. 

PACS nos.: 25.20.Lj , 12.39.Fe 



1 Introduction 

One of the major goals in nuclear physics is the understanding of isospin symmetry vi- 
olation related to the light quark mass difference rriu — rrtd and virtual photon effects. 
Although the light quark mass ratio deviates strongly from unity, nid/mu ~ 2 0] and 
one thus could expect sizeable isospin violation, such effects are effectively masked since 
iT^d — TTT-u -C A, with A the scale of the strong interactions (which can be chosen to be 
47rF7r or 1 GeV or the mass of the p). To assess the isospin violation through quark mass 
differences, precise measurements and accurate calculations are mandatory. As pointed 
out by Weinberg a long time ago [0], systems involving nucleons can exhibit such effects 
to leading order in contrast to the suppression in purely pionic processes due to G-parity. 
Pion-nucleon scattering or neutral pion photoproduction are best suited for such inves- 
tigations since a large body of (precise) data exists for various isospin channels. While 
there is still some discrepancy about the low-energy tt~p elastic scattering data, very 
precise neutral pion photoproduction experiments have been performed over the last few 
years at MAMI and SAL [Q. Together with the new data on 7p — >• Tc'^n, a sensitive 
test of isospin symmetry could be performed if one would have information on the ele- 
mentary neutron amplitude 'yn -^ ir^n. This amplitude can only be inferred indirectly 
from reactions involving few-nucleon systems like the deuteron or ^He. In this paper, we 
concentrate on coherent neutral pion production off deuterium in the threshold region 
and study the sensitivity of the deuteron electric dipole amplitude Ed to the elementary 
neutron amplitude, -Eo+". The framework to do this is baryon chiral perturbation theory, 
which is reviewed in great detail in Ref. 0. It is the effective field theory of the stan- 
dard model at energies below 1 GeV and allows one to explore in a systematic fashion 
the strictures of the spontaneously and explicitly broken chiral symmetry underlying the 
fundamental color gauge theory of the strong interactions, QCD. While originally formu- 
lated for the Goldstone boson sector, the machinery can be extended straightforwardly to 
include single-baryon processes, treating the baryons as very heavy, static sources [Q . 
In few-nucleon systems, a complication arises due to the existence of shallow nuclear 
bound states and the related infrared singularities in Feynman diagrams evaluated in 
the static approximation. One way to overcome this is to adapt the rules and use chiral 
perturbation theory to calculate an effective potential, which consists of the sum of all 
A-nucleon irreducible graphs 0. The S- matrix, which of course includes all reducible con- 
tributions, is then obtained through iteration by solving a Lippmann-Schwinger equation. 
Several generic features of nuclear physics, in particular the relative strengths of two- and 
few-body forces [Q ||TD[ |TT| or the dominance of soft-pion exchange currents |T2|, have been 



shown to arise naturally in this approach. For calculating scattering processes involving 



light nuclei, Weinberg |jT3[ proposed to use chiral perturbation theory to generate the ir- 
reducible kernel and combine these with very accurate phenomenological external nuclear 
wave functions. Although nuclear wave functions are calculable in chiral perturbation 



theory |Ty], they have not yet been determined with the accuracy matching that of the 



phenomenological models. Therefore, we use the input from the successful phenomenolog- 
ical boson-exchange potentials which differ mostly in their short-range interaction parts. 
Using a large variety of these potentials allows one to assess to which degree of accuracy 
one is sensitive to the chiral symmetry constraints used in determining the irreducible 
scattering kernel. 

In this paper, we consider pion photoproduction on the deuteron at threshold extending 
the previous work by Beane et al.|T^. We calculate the invariant threshold amplitude 



to fourth order in small momenta and/or pion masses. This is mandated by the fact 
that only to this order in the chiral expansion one can achieve a satisfactory description 
of the elementary process jp -^ n^p [|^[|1^. As already stressed in Ref. ||14||, to third 



order in small momenta, the interaction kernel for neutral pion photoproduction off the 
deuteron has no undetermined parameters. Unfortunately, although the same is true in 
the single nucleon sector, evidently the amplitude there converges slowly at best and 
thus the single-scattering contribution was treated as phenomenological input. Since 
that paper was published, the situation has considerably improved and we are now in the 
position to consistently include the predictions for the single-nucleon sector. In particular, 
the precise 0{q^) (where g is a generic symbol for a small momentum or meson mass) 
chiral perturbation theory calculation [T^ [|l^ predicts a very large neutron electric dipole 
amplitude, larger in magnitude than the corresponding proton one. As we will show, to 
fourth order in small momenta the interaction kernel can still be given entirely in terms 
of known parameters, and thus we will be able to calculate the sensitivity of the threshold 
cross section to the elementary neutron-vr" amplitude. We mention that charged pion 
photoproduction at threshold has been calculated accurately (to four orders) in chiral 
perturbation theory |l^ . 



Since neutral pion photoproduction has already been studied extensively in more con- 
ventional nuclear approaches based e.g. on meson-exchange models (see e.g. refs.|l^- 



), it is worth emphasizing why we use chiral perturbation theory here. As first stressed 
by Weinberg |]l3l, chiral perturbation theory allows one to systematically construct the 
many-body interactions between nucleons, pions and photons following the power count- 
ing rules. For example, in the case of elastic pion-deuteron scattering the phenomeno- 
logically dominant rescattering contribution was shown to be the leading three-body 
interaction in the chiral expansion thus leading to deeper theoretical insight. Here, the 
phenomenological approaches to the single neutral pion production amplitude off pro- 
tons based on Born terms are at variance with the data by many standard deviations 
quite in contrast to chiral perturbation theory. To tackle the problem of how sensitively 
the deuteron threshold cross section depends on the elementary neutron-7r° amplitude can 
therefore only be addressed in a framework which allows one to systematically include and 
order the various contributions arising from single and multiple scattering processes. It 
goes without saying that this approach is only useful for the threshold region. We remark 



that contributions from the A-isobar are encoded in some of the low-energy constants 
appearing in the single scattering matrix elements. 

This paper is organized as follows. In section 2 we briefly review the effective La- 
grangian underlying the calculation and the standard power counting formulas. In sec- 
tion 3 the calculation of the various contributions to the threshold S-wave amplitude is 
outlined. Section 4 contains the results and discussions thereof. The appendices include 
our conventions and give some more details on the calculations. 

2 Effective field theory 

In this section, we briefly discuss the effective chiral Lagrangian underlying our calcula- 
tions and the corresponding power counting. For more details we refer to the review 
and to the paper by Beane et al. |T^ . 



At low energies, the relevant degrees of freedom are hadrons, in particular the Goldstone 
bosons linked to the spontaneous symmetry violation. We consider here the two flavor 
case and thus deal with the triplet of pions, collected in the matrix U{x) = C,'^{x). It 
is straightforward to build an effective Lagrangian to describe their interactions, called 
Ct^tj. This Lagrangian admits a dual expansion in small (external) momenta and quark 
(meson) masses as detailed below. Matter fields such as nucleons can also be included 
in the effective field theory based on the familiar notions of non-linearly realized chiral 
symmetry. These pertinent effective Lagrangian splits into two parts, C-j.n' and Cj^n, 
with the first (second) one consisting of terms with exactly one (two) nucleon(s) in the 
initial and the final state. Terms with more nucleon fields are of no relevance to our 
calculation. The pertinent contributions to neutral pion photoproduction at threshold 
are organized according to the standard power counting rules, which for a generic matrix 
element involving the interaction of any number of pions and nucleons can then be written 
in the form 

M = g'^^(g//i), (1) 

where /z is a renormalization scale, and i/ is a counting index, i.e. the chiral dimension 
of any Feynman graph, u is, of course, intimately connected to the chiral dimension di 
which orders the various terms in the underlying effective Lagrangian (for details, see 0). 
For processes with the same number of nucleon lines in the initial and final state {A) , one 
finds§ 

z/ = A-A-2C + 2L + J2yi^i 

i 

A, = di + ni/2 - 2. (2) 

where L is the number of loops, Vi is the number of vertices of type i, di is the number of 
derivatives or powers of M^r which contribute to an interaction of type i with Ui nucleon 



fields, and C is tlie number of separately connected pieces. This formula is important 
because chiral symmetry places a lower bound: Aj > 0. Hence the leading irreducible 
graphs are tree graphs (L = 0) with the maximum number C of separately connected 
pieces, constructed from vertices with Aj = 0. In the presence of an electromagnetic field, 
this formula is slightly modified. Photons couple via the electromagnetic field strength 
tensor and by minimal substitution. This has the simple effect of modifying the lower 
bound on A^ to Aj > —1, and of introducing an expansion in the electromagnetic coupling, 
e. Throughout, we work to first order in e, with one exception to be discussed below. In 
what follows, we will work within the one-loop approximation to order q^ (notice that we 
refer here to the chiral dimension used to organize the various terms in the calculation 
of the single-nucleon photoproduction amplitudes), which extends the 0{q^) calculation 
of ref. [l^ . Such a higher order calculation is mandated by the fact that to order q^ the 
single-nucleon neutral pion photoproduction amplitudes are too inaccurate. Furthermore, 
one would like to see how big the corrections to the leading order three-body interactions 



calculated in ref.|T^ are. In terms of the counting index u, in ref. [|T^] all terms with 
z/ = 4 — 3A = —2 and z/ = 5 — 3A = —1 were included. We go one order further, i.e. 
we add all terms with u = 6 — 3A = (remember that A = 2 and we have exactly one 
photon coupling with Aj = —1). Consequently, the effective Lagrangian consists of the 
following pieces: 

r - r(^) -L r(i) -l r^^) , A^) ^(4) ^(o) ^(i) /o\ 

where the index (i) gives the chiral dimension di (number of derivative and/or meson mass 
insertions). The form of C^} + C\^l, is standard. The terms from £^^ + £^jy contributing 
to the single-nucleon photoproduction amplitudes are given in ref. |]T3|. Concerning the 
terms from £^j(f and C^i^., we will discuss the pertinent ones below.Q With these at 
hand, we have to calculate tree and loop graphs for the single-scattering amplitude and 
tree graphs with exactly one insertion from C)^l^ for the three-body interactions between 
pions, nucleons and photons as well as tree graphs with one insertion from C}/j^. As 
we will show later, the only new coupling constants appearing at 0{q'^) are related to 
the single-nucleon sector and have already been determined in refs.|jl5||[T^. After these 
general remarks, let us now turn to the actual calculations. 

3 Anatomy of the calculation 

In this section, we outline how the various contributions to the electric dipole amplitude 
Ed are calculated. The pertinent details are relegated to the appendices, including our 
conventions and definitions of the S-matrix. 



^Of course, terms from £^^ also appear in the single nucleon calculation of refs.jIq][E6J and from C\/j^ 
in the nuclear force calculation of ref. lllli . 



3.1 Threshold cross section 

Consider the reaction 'y{k) + d{pi) -^ vr°(g) +d{p2) in the threshold region, g ^ 0. For real 
photons [e- k = 0, with e the photon polarization vector) the S-wave photoproduction 
amplitude can be expressed in terms of two multipoles as, 

Md = 16z7r(md + M^) J d^p(l)}{p) |a^i -(^i + ^2) ■ e 

+ M2^[t^i-ka2-iexk) + il^2)]'j(j),{p-k/2) . (4) 

The differential cross section at threshold takes the form (see appendix A) 



|fc| da 
\q\ dQ 



= \Ei . (5) 

|q|=0 O 



In what follows, we will discuss the chiral expansion of the dipole amplitude E^. In 
appendix A we discuss the relation to the commonly used amplitudes treating the deuteron 
as an elementary particle. 

3.2 Single scattering contribution 

The single scattering contribution is given by all diagrams where the photon and the pion 
are absorbed and emitted, respectively, from one nucleon with the second nucleon acting 
as a spectator (the so-called impulse approximation). One only has a contribution to the 
A4 1 amplitude of the form 

_ ^k. fd'pp\ {Pf^ + Pf-) ^}ip) e- J0,(p - k/2)] , (6) 

m J 2 ■' J 

evaluated at the threshold value 

\k\ = fcthr = M^o - -^ = 130.1 MeV , (7) 

2md 

and with J = {ai + 02)12. A number of remarks concerning Eq.(|^) are in order. It 
is important to differentiate between the -n^d and the 7r°iV (A^ = p, n) center-of-mass 
(COM) systems. At threshold in the former, the pion is produced at rest, it has, however, 
a small three-momentum in the latter [|1^. Consequently, one has a single-nucleon P- 
wave contribution proportional to the elementary amplitudes P^ ^ and PI " as defined in 
|TB|. We use the P-wave low-energy theorems found in that paper. 



Pl P = 0.480 |g I GeV"^ Pf " = 0.344 |g | GeV" 



with 

g = /i(l-/i)p-/i2m(l-5/i/4)A;/2 , (9) 

with /x = MT^/m and p is the nucleoli three-momentum in the vr^d COM system. The 
derivation of Eq.@ is sketched in appendix B. We have checked that dropping the terms of 
order M^ in Eq.(^ does not alter the results within the precision given. Furthermore, we 
neglect the energy dependence of the elementary 7r°p and vr'^n S- and P-wave amplitudes 
since the pion energy changes only by 0.4% for typical average nucleon momenta in the 
deuteron (see appendix B). We take for the elementary S-wave pion production amplitudes 
the predictions from the 0{q^) chiral perturbation theory calculation ||16|| , 



^0 



^J/ = -1.16 X 10"7M^+ , EJ/ = +2.13 X 10"7M^+ . (10) 

In that calculation, the dominant isospin breaking effect, namely the charged to neutral 
pion mass difference, which is almost entirely of electromagnetic origin, has been taken 
into account. A fully consistent calculation including all effects from virtual photons 
and the quark mass differences is not yet available. The result for Eqj^ is based on the 
assumption that the counter terms entering at order q^ are the same for the proton and 
the neutron apart from trivial isospin factors (for a detailed discussion, see ref.[^). We 
also note that only the first process has been measured, the recent determinations from 
MAMI and SAL give 

Ef/ = (-1.31 ± 0.08) X 10-VA'4+[3] , Eq? = (-1-32 ± 0.08) x lO^V^^+l • (H) 

Furthermore, the P-wave LET for neutral pion production off protons has been shown to 
hold within 3% in Ref.[|]. Using the Argonne V18 H, the Reid soft core (RSC) |23|, the 



Nijmegen p4[ and the Paris |^ potential, we evaluate E^^ and find 

Ef = 0.36 X 10^7^^+, (12) 

with an uncertainty of 5E^^ = 0.05 x 10~'^/Mn-+ due to the various potentials used. The 
P-wave contribution amounts to a 3% correction to the one from the S-wave, i.e. it 
amounts to a minor correction. The sensitivity of the single-scattering contribution E'J* 
to the elementary neutron-7r° amplitude is given by 

^^« = 0.36 -0.38- (2.13 - EJ+") , (13) 

all in units of 10~^/Mn-+. Consequently, for £'0+" = 0, we have E^^ = —0.45 which is of 
opposite sign to the value based on the chiral perturbation theory prediction for -Eo+"- 
If one were to use the empirical value for the proton amplitude, the single-scattering 
contribution would be somewhat reduced. 



Finally, we remark that the single-scattering contribution given in ref.|Tj] is very dif- 
ferent. This is due to the following. First, in that paper a factorized form for the single- 
scattering contribution was used. 



ET = 1+mJZ ^ ^^°^' + ^'"-^^ S,{h^j2) = 0.41 X 10-VM., 



(14) 



where S'd(A;thr/2) is the deuteron form- factor, 

Sdih^r/2) = f d'p (P}ip) (Piip- 4hr/2) = 0.79 . (15) 



Note that the isospin factor 1/2 was inadvertently omitted in [|I4|. Such a form is strictly 



correct only for the S-wave part of the deuteron wave function. Switching off the D-wave 



component of the deuteron wave function, we get the same result as in jT^. The main 



reason for the large numerical difference of our single-scattering contribution to theirs 
can be traced back to the use of the then accepted empirical value of -Eo+^ = (— 2.1 ± 
0.2) X 10^^/M^+ and the one from the incomplete low-energy theorem for 'yn -^ vr^n, 
EJ+" ~ 0.5 X 10-VM^+ in 0] . 

3.3 Three— body contributions at order q^ (i^ < —1) 

Although the 0{q^) (corresponding to the counting index u = —2 and —1) three-body 



contributions (meson exchange currents) have been worked out in ref. [0, we will give 
them here for completeness. In the Coulomb gauge {e ■ v = 0, with v^ the nucleons' 
four-velocity) and at threshold, only the two diagrams shown in Fig. 1 contribute. In 
momentum space, the corresponding structures are 



]_ ^. (g - k)^ qj 



4, b-- ...'v::'':l.. , m 



where q = p — p' with p and p' the nucleon cms three-momenta in the initial and the final 
state, respectively. These expressions can be Fourier transformed into coordinate space 
easily and one finds that they give a contribution to the J^i amplitude defined in Eq.(^ 
]I^ . Evaluating these with standard deuteron wavefunctions gives. 



where the indices S, SD, D refer to the contribution from the S-wave, the mixed and the 
D-wave part of the deuteron wavefunction. We find for the three-body (tb) contribution 
at order q^ 

^tb,3 ^ (_i 90^ -1.88, -1.85, -1.88) x 10"V^^+ , (18) 

for the Argonne V18, RSC, super soft core (SSC) |2^ and Bonn potential ||27 



m 



order, and using g^ = 1.32 as determined from the Goldberger-Treiman relation (for 

8 



Qt^isi = 13.4) (to be consistent with the calculation of the elementary amplitudes [0). 
Note that we use the full (energy-dependent) model of ref. [^ throughout. The number 



for the Bonn potential differs from the one given in ref.|Tj] for three reasons. First, the 
overall sign of the contribution from graph b was incorrectly given and second, we find 
that the SD- and D-wave contributions from this graph substantially alter the pure S- 
wave contribution. Also, the D-wave contributions to graphs a and b were incorrectly 
evaluated (the corrected formulae are given in appendix C). If we, however, only retain 



the S-wave part of the deuteron wavefunction, we recover the result of Eq.(26) of ref.|T4 
(modulo the sign and within the numerical precision). 

3.4 Three— body contributions at order g^ (v = 0) 

At this order, we have to consider tree graphs with exactly one insertion from the (chiral) 
dimension two pion-nucleon Lagrangian, which has the general form (in the isospin limit 

rriu = rud) 

-^[-^^ S:] [(1 + «:.)/+ + li^s - «:.)Tr(/+ )] (19) 

+ci Tr (x+) + (c2 -^yvuf + csu-u+ (^c^ + ^) [^^ S^ju^u, ^H, 

where m is the nucleon mass and /^ = e{^'^Q^+^Q^^)F^^. F^^ is the electromagnetic field 
strength tensor and 5''^ is the covariant spin-operator. H denotes the large component of 
the nucleon spinor (for more details, see ref.p). The terms in £^^ fall into two classes, 
the first one with fixed coefficients due to Lorentz invariance[0 and the second one with 
low-energy constants Ci, . . . , C4 and two other constants that are directly related to the 
isoscalar and isovector nucleon anomalous magnetic moments, k^ = Hp — Hn, ^^s = i^p + i^n- 
The pertinent three-body interaction diagrams are shown in Fig. 2, the blob characterizing 
an insertion from /^jtTV- I^ principle, all the terms appearing in Eq.(|l9l) can contribute. 
In addition, there are the recoil corrections to the graphs in Fig. 1. 

Consider first insertions proportional to the low-energy constants Ci^2,3,4- One can show 
that all of corresponding terms vanish at threshold either due to the selection rule S -q = 
or due to some isospin factor of the type 6°'^ e"'^^. We are thus left with insertions from the 
terms ~ l/2m, ~ (7^/2m and ~ Ky^g. All of these can be classified in momentum space 
according to the following structures. 



1 qi qi qj 
^2 ' ^2 ' ^2 ' 
q q q 

iq-k)i qj 


qi{p + p')j {q-k)i{p + p')j 

{q-k)iqj{p'^ -p'^) {q-k)i 


(q - k); qj 


(g - k)i {q - k)j 


\q-ky + Mlf/^ 


' [(g-'_^)2 + M2]3/2 ■ 



(20) 



All these operators can straightforwardly be transformed into coordinate space and the 
angular integrations can be performed analytically leading to the form given in Eq.(P and 
one is left with one simple radial integration. Only the second operator in the second line 
is more easily evaluated in momentum space. Whenever possible, we have performed both 
types of integration as a check on our numerics (for details, see appendix C). In fact, there 
is a technical problem with the contribution from graphs b, 1 and p. Counting powers 
of momenta, one sees that they are divergent. This is due to the fact that in the chiral 
expansion one has truncated the pion-nucleon form factor and thus does not suppress 
the high-momentum components. This phenomenon also occurs in the calculation of the 
NN potential in chiral perturbation theory []TU|. As it was done there, we introduce an 



additional Gaussian cut-off factor of the form 

F{f) = exp{-f/Al} , (21) 

where the cut-off A^ varies between the mass of the p and 47rF^ = 1.2 GeV (since it can be 
related to the mass scale of the heavy particles that are integrated out from the effective 
theory). Notice that the divergencies appearing in other graphs all cancel and thus no 
further regularization is needed for those. For energy-independent NN potentials as used 
here, the contribution from the graphs n + r should be omitted for the reasons explained 



in detail in ref.[|lj]. We will, however, also calculate their contribution to get a very rough 
estimate of some of the uncertainties due to higher order graphs. Furthermore, we remark 
that the kinematical 1/m corrections to graph la vanish and the ones for graph lb are 
taken care of by shifting the value of fcthr as given in Eq.(|^). Finally, appendix D contains 
some details on calculating the time-ordered diagrams in the heavy baryon approach. 
We find for the three-body (tb) contribution at order q^ 

Ef^ = (-0.25 , -0.23 , -0.27) x 10-^/M^+ , (22) 

for the V18, Reid and SSC potentials, in order, and setting A^ = 1 GeV. These amount 
to corrections of the order of 15% to the q^ three-body terms. If one varies the cut- 
off Att from 0.65 to 1.5 GeV, the three-body contribution using the V18 potential varies 
between —0.24 and —0.29 (in canonical units), which is a modest cut-off dependence. The 
result Eg . (^2]) is comforting since it shows that the chiral expansion of the three-body 
contributions is well under control. Finally, we remark that if one were to include the 
contribution from the recoil graphs n + r, the numbers given in Eg . (|2^) would change by 
less than 5 per mille. 



10 



3.5 Contribution from four— nucleon operators (z/ = 0) 

In the O(g^) calculation of ref. [Q, no four-nucleon operators contributed to the deuteron 
electric dipole amplitude. Using Eq.(|D, we see that for z/ = we can have exactly 
one insertion from jC^vtv- I^ F^S- 3a,b, we show the two graphs which in principle can 
contribute to pion photoproduction at threshold. Clearly, diagram 3a is only relevant for 
charged pion photoproduction. Graph 3b also vanishes for neutral pions, as the following 
argument shows. This diagram stems from the one shown in Fig. 3c by applying minimal 
substitution. The latter one contains the pion covariant derivative, V^ = d^ji jF^, + . . . 
which upon minimal substitution takes the form 

diJ^a -* diJ^a -ieA^ Qab T^b , (23) 

with 

Qab = i Gabs , (24) 

the pion charge matrix and A^ the photon field. Consequently, d^n^ — ^ d^n^ and this 
type of term can only contribute to charged pion photoproduction. We therefore conclude 
that for threshold neutral pion photoproduction to order g^ (or up to counting index 
z/ = 0) there is no contribution from any four-nucleon operator and thus no new, a priori 
undetermined coupling constants appear. In the case of charged pion photoproduction, 
this would be different. 

4 Results and discussion 

Since as we have shown in the previous sections, neither the single scattering nor the 
three-body corrections depend on the potential chosen, we will present here results based 
on the Argonne V18 potential. The chiral expansion of the electric dipole amplitude Ed 
takes the form 

Ed = E^' + Ef^ + Ef^ 

= (0.36 - 1.90 - 0.25) X 10^7^4+ = (-1-8 ± 0.2) x 10"V^^+ • (25) 



It is difficult to estimate the theoretical uncertainty, the value given in Eq.(|25D being 

— 

an educated guess, obtained as follows. We allow Eq^^ to vary between —1. and —1.5, 
so that for a fixed value of -Eo+" ~ 2.13, E^'^ varies between 0.2 and 0.4 (all numbers 
in canonical units). The results for the three-body contributions are stable at order q^, 
where as we assign a conservative uncertainty of ±0.1 to the corrections at next order 
due to the cut-off dependence. Clearly, this estimate does not include contributions from 
higher order effects, which might give rise to a larger uncertainty. Our final value is 
considerably smaller than the g^ estimate reported in fl^ for the reasons given above. To 
see the sensitivity to the elementary neutron-7r° amplitude, we set the latter to zero and 

11 



find Ed = —2.6 X 10~^/M^+ (V18 potential) which is considerably different from the chiral 
perturbation theory prediction, Eq.(^), i.e. the S-wave cross section would differ by a 
factor of two. For other values of -Eon-"' ^d can be calculated from Eg. ([T3|) . Obviously, 
the sensitivity to the neutron amplitude is sizeable and is not completely masked by the 
larger charge-exchange amplitude as it is often stated. 

On the experimental side, neutral pion photoproduction off deuterium was studied 
by a group at Saclay ^9| and later reanalyzed in ref.[Q. In these papers, the S-wave 
amplitude S2 was extracted. We remark that this analysis relies heavily on the input 
from the elementary n^p amplitude (to determine an unknown normalization factor) and 
therefore should only be considered indicative. The S- and P-wave multipoles used in 
PPI are roughly consistent with the new determinations from Mainz and Saskatoon. The 
amphtude S2 is related to our Ed via 



Using Sd = 0.79 and the value of £2 = (-4.1 ± 0.4) x 10-^/M^+ |0l, we have as the 
"empirical" value 

^f P = (-1.7 ± 0.2) X 10'yM^+ , (27) 

taking the same sign as given for 82- Since we could not trace back what the exact value 
of 5*^ used in the Saclay analysis was, we took the value based on the modern potentials 
evaluated above. Notice that in the older papers of the Saclay group, a larger value of £2 
was obtained, based on different input for the proton amplitude (which is at variance with 
the new data from Mainz and Saskatoon). This was used to deduce the empirical number 
quoted in ^^. Obviously, this experiment should be redone at a 100% duty cycle tagged 
photon facility, with the same accuracy as was done for the process 7p -^ n^p at MAMI 
and SAL and with a refined theoretical analysis as is available now. The empirical number 
Eq.(^) agrees nicely with the theoretical one, Eq.(^). However, we remind the reader 
about all the assumptions going into the extraction of this "experimental value" . Clearly, 
only a more precise experimental determination of Ed can tell whether this agreement is 
of significance. On the other hand, if one would find a discrepancy, one would either have 
to reassess the calculation of the elementary amplitudes by including dynamical isospin- 
breaking effects and/or study in more detail the wave function dependence of the order q^ 
three-body corrections. Finally, two more remarks are in order. First, the sensitivity to 
the neutron amplitude has recently been studied in more conventional approaches [^ |^ . 



Second, an experiment has been approved at the Mainz Microtron [^ to measure the 
threshold cross section for coherent neutral pion electroproduction off deuterium at a 
photon virtuality of k"^ = —0.075 GeV^. 
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A Invariant amplitudes for the deuteron 

If one considers the deuteron as an elementary particle, the S-wave photoproduction 
amplitude can be written as 

Md = 8TT{md + M^)2te-JEd + 0{q) , (A.l) 

with rriii the deuteron mass and J = L + S the deuteron total angular momentum, not 
simply the sum of the proton and the neutron spin operators. This is the form commonly 



used in analysing the data. Following the conventions used in |]IJ], the slope of the 
differential cross section at threshold takes the form 



\k\ da 

1^ 'dn 



\M, 



.1=0 647r2 Urnl + Ml, + |fc|) (m, + M^.] 



1 |A<dP 



647r2 {md + M^of ' 
Summing over the final and averaging over the initial states leads to 

2 

191=0 2 J -, ^ MM\X 



(A.2) 



\k[ da_ 

J^ dn 



- Y. \<JM'\ex-J\JM> Ed 

^ M Ml 



|2 



\{2J+l)\Ed\' = l\Ed\'' (A.3) 



2J+1 3 ' '' ' 3 



where M, M' are magnetic quantum numbers and A the helicity index of the photon. We 
have made use of the Wigner-Eckhardt theorem and the fact that < J|| J || J >= \/2 for 
the deuteron. 

B Two-body to three-body center-of-mass 

In this appendix we sketch the derivation of the transformation from the '-j-d center-of- 
mass (COM) system to the 7-iV COM. We are interested in the kinematics of the process 
'yNiN2 —>■ ttNiN2, where the nucleons, Ni and A'"2, are sewn to the deuteron wavefunc- 
tions. Our 3-body corrections are evaluated in the j-d COM whereas the single scattering 
corrections which take into account the scattering of the photon on the individual nu- 
cleons have been calculated in the 7-A^ COM. It is therefore necessary to construct the 
Lorentz transformation which boosts the single-scattering corrections to the 'j-d COM. 

We label energy-momenta as pi, p2 and k for A^^i, N2 and 7 in the initial state, respec- 
tively, and p[, p'2 and q for Ni, N2 and vr in the final state, respectively. The energy- 
momentum of the deuteron is given by p^ = pi + P2 and p'^ = p\ -|- p'2 in the initial and 
final state, respectively. The Fermi momentum, p, can then be defined by p\ = p^/2 + p 
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and p2 = Pv/2 ~ P- Note that with this labelhng the initial deuteron wavefunction has 
the argument (pi — p2)/2=p. 

The 7-(i COM is defined by p^ + /c = and the 7-iV2 COM is defined by p2 + ^ = 0. 
The velocity of the 7-A'^2 system in the 7-(i COM is 

P=^ (B,l) 

where i?27 = E2 + k. An arbitrary four- vector, {E,p), in the '-f-d COM is expressed in the 
7-A''2 COM by the matrix equation 





The *-superscript indicates the 7-A^2 COM and p = p^ + P\\P- We also have pl_ = p±. 
We can now use this transformation to find k and q in both frames. In the text we use a 
convention in which the initial deuteron wavefunction has argument p — k/2. Here this is 
achieved by the substitution p — > —p + k/2, such that p2 = p — k and pi = —p, since N2 
must carry all the photon momentum. At the threshold point we find 

k^ = ko — k • p/m , k* = k — (kQ/m) p , 

E: = E^[1- f/(2m')] , q* = -E^p/m , (B.3) 

where we have included the first non- vanishing 1/m corrections. Note that the pion energy 
is corrected only at order l/rn^. Calculating these corrections with an average nucleon 
momentum in the deuteron obtained by means of the uncertainty principle, (p ) = 46 MeV, 
leads to the 0.4 % shift mentioned in the main text. Of course, such an crude estimate 
does not properly account for binding energy effects and should be sharpened eventually. 
In the 7-A^2 COM we have the multipole decomposition 

T-e = ia- e{Eo+ + k* ■ q*Pi) + ia ■ k* e- q* P2 + (g* x A;*) ■ ePs- (B.4) 



In terms of the transformed variables it is straightforward to find that the S-wave multi- 
pole is modified by 

Eo+^Eo+--^pA (B.5) 

m 

where we have left out terms that vanish upon integration over the Fermi momentum, 

and it is understood that the argument of Pi is a transformed pion momentum. 
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C Nuclear matrix elements at orders (^ and q' 



In this appendix, we give the coordinate space representations corresponding to the di- 
agrams shown in Fig. 2. Whenever appropriate, we combine graphs. Graph d vanishes 
because of isospin and e and f are proportional to v ■ {q — k) which is zero at threshold to 
the order we are working. 



a+k+o: 



b+l+p: 






(c.i: 






+ - / rfV(CTi+CT2) -f 

TT 






-2 0,^^(^ + 2 0/^(^^^) 









r 



/2 



t-r—i d €■ (r + r 



(r*-!- r ') . f d (r*+ r ') ■ k 
r dr 2 



ik-1 



[r 



c+i+j: 



— ?- 



5egA 



57r 



(! + /€,) d'r(l>*{r)03 



Yi{M^r) 



3.^*, ^^ VllM^O^.^f 



e 2 r 



(C.2) 



(C.3) 



r2(M^r) Fi(M.r) 8 



g+h: 



m+q: 



+ / dM*{r) e-'^{ -Oi ^^^5^ + O, ( '^"7"^ - 2 



' rf>0*(r)e ' 2 J_Oi — + 0: 



(ViO /'j/Y' 



AtxF^ 



^{/<;Vr(r)e-* 



-0.^1(^.0, 



Y^iM^r) Yi{M^r) d 



(C.4) 



"^ C/T" 



dV^-Ml^o.^ii^ + o,^^*^^-*^,,) ^ . (c,5) 



n+r: 



with 



SeglM^ 



rf^r0*(r)e*— O4 

Yi{M^r) : 
F2(M.r) : 



dKo{M^r)_^^d^Ko{M^r) 



dr 

M^ + - I e 
r 



Q^2 



(C.6) 



-M^r 



, .2 3M^ 3 

^' + V^ + ;:^ 

°° z sin 2; 



-M^r 



(^2 + M2r2)3/2 ' 



(C.7) 
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and 



O4 



0-1 + 0^2) ■ e*, 
f?i + f?2) ■ re- f , 

(Ti ■ f?- (f X (T2) + (1 ^ 2) ) , 



(CJ 



Performing the angular integrations, all the spin-dependence can be expressed in terms 
of the spin vectors a\^2 and the photon three-momentum k as given in Eq.(|^). The second 
operator in Eq.(^ stems in part from O2, O3 and O4. The deuteron wave function 0(r ) 
is given by 

,(,,) 1(^,1.,,,) EM). ,c.9) 

V47r V r V8 r ) 

in terms of the S- and D-wave functions Uir) and W{t\ respectively, and 

S\2 = 3(ai ■ f){a2 ■ f) - (fi ■ (72 (C.IO) 

is the tensor operator. The wave function is normalized to one, /q°° dr (f/^ + W"^) = 1. 

For completeness, we also give the corrected version of the g^ contributions correspond- 
ing to graphs la and lb. The angular integrations are performed as outlined in ref.||I^ 
using the corrected relation 

SI, ( J • e)Su = 4 [3{e- f){J-f)-2 (J ■ e)] . (C.ll) 



This gives: 
a: 



iegAMj^m r dr 



vrFi 



/dv ( 1 
— - \ U'^ 0? sin a H ^ f/ ly (3 sin a — 3a cos a — a^ sin a) 



+ -W'^ (3 sin a — 3a cos a — 2a^ sin a] 



(C.12) 



b: 



iegAMT^m 



ttF^ 



dz dr- , „ 
J rb-^ 



1 



t/^ b^ smb + ^UW{3 sin 6 - 36 cos 6 - 6^ sin 6) 



v^ 



+ -Vr2 (3 sin 6 - 36 cos 6 - 26^ sin 6)] 
dz (ir YifmV)— ■ (sin6 — 6cos6) ( f/ H ^ly 

7o 7 ^ ^ 63 ^ ^ V ^/2 



(C.13) 



3|A; 



dz iz — 11 



-m r , 1 

rfr—— f?7iy + ^iy2) ((3 -6^) sin 6 -36 cos 6) 



with 



kr , , / 1' 
a = y, 6 = fcr(^z--,, 



m' = M^ a/z (2 - z) . 



(C.14) 
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D Time orderings in heavy fermion formalism 

In this appendix we establish a simple method for extracting irreducible time ordered 
graphs in baryon chiral perturbation theory. Consider the reducible relativistic Feynman 
graph shown in Fig. 4. Since there are 3 interaction vertices we expect 3! time order- 
ings. The time ordered decomposition of Fig. 4 is displayed in Fig. 5. There are several 
important points to note: 

• Graphs (4), (5) and (6) are 1/m corrections and so the sum of (1), (2) and (3) correspond 
to the Feynman graph of Fig. 4 evaluated in the heavy fermion formalism (HFF). 

• Graphs (2) and (3) are reducible and therefore graph (1) is the irreducible subgraph 
that we are after. 

The propagator structure of Fig. 4 is 
i i{;^ + m) i 



k'^ — Ml p2 _ ^2 2ci; V /co — c^ A^o + ^ 



% f'-fQEn-l-p + m 7o^n + 7-p-m ^. 



2En V Po — En Po + En 



where En = \Jp^ + m^ and uj = yk? + M^. We are interested in positive frequency pions 
and nucleons and so 

AujEn \Ko — UJ Pq — En / 

This expression can be partial fractioned to give 

-7oE„ + 7-p-m 1 / 1 1 \ .j^ gx 

AujEn ko + Po - En - LOKko - UJ Pq- En)' 

In order to show that the two pieces inside the parentheses are in correspondence with 
(1) and (2), we evaluate (1) and (2) using old-fashioned time ordered perturbation theory. 
Here we attach a photon and a pion line to our generic interaction vertex and focus on 
the distinct time slices, labelled {i) and {ii) (see Fig. 6). The rules are simple: each 
distinct time slice — with respect to the interaction vertices — corresponds to a single 
energy propagator. We find 

(1) oc ^ 



[(^1 + E2 + E^) -{E[+u + En + E^)] 

1 
[(^1 + E2 + E^) -{E[ + uo + E2 + E^)] 
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[ii) (D.4) 



"" [{E^ +E, + E,) - {E, + E„ + E^)] ^''^- ^^-^^ 

In relativistic notation we have ko = Ei — E[ and pq = E2 — E^ — E.^. (Note that energy 
is not conserved at the vertices in the time ordered formahsm.) It then follows that 

1 1 

(1) oc 1 F X 



ko + Po - En - uj ko-u 
(2) oc ^- X ^-—. (D.6) 

ko+po- En-UJ Po- En 

We are now in position to establish a rule for extracting the irreducible time ordered 
subgraph, (1), in HFF. The propagator structure of Fig. 4 in HFF is 

(D.7) 



k^-M^ vk 

where p = mv + k. In the static limit v ■ k —^ kQ=0, which renders (D.7) singular. In this 
limit En = m and po = m — ko. Therefore, decomposing the propagators gives 



i ^ { ^ ~^ 



X 



k'^ — Ml ko 2uj\ko — to ko + to/ po — En 

^ ^ (-)T.O. 



2uj {ko - uj){po - En 



luj ko \ Pq — En — ijj\ko — ijJ Po — E. 
1 1 1 



reducible (+) T.0. + (-) T.O. 



2uj {ko +po- En-uj) {ko - to) 

(D.8) 

where (±) refers to the frequency of the time ordering and (D.6) was used in the last line. 
The first term is graph (1) of Fig. 5, the irreducible (+) T.O. piece. It is, of course, well 
behaved in the static limit. In the limit /cq ^ we have 

1 1 
irreducible {+) T.O. = — - = tj:^ (D.9) 

2^ 2(F + M^f 
and so the irreducible time ordering in a HFF Feynman graph of type Fig. 4 can be 
extracted by making the following replacement in the HFF propagators: 

(D.IO) 



k^-M^vk 2{k- + Mlf'' 



19 



References 

[1] J. Gasser and H. Leutwyler, Nucl. Phys. B94, 269 (1975). 

[2] S. Weinberg, Trans. N.Y. Acad. Sci. 38, 185 (1977). 

[3] M. Fuchs et al., Phys. Lett. B368, 20 (1996). 

[4] J. Bergstrom et al., Phys. Rev. C53, R1052 (1996). 

[5] V. Bernard, N. Kaiser, and Ulf-G. MeiBner, Int. J. Mod. Phys. E4, 193 (1995). 

[6] E. Jenkins and A.V. Manohar, Phys. Lett. B255, 558 (1991). 

[7] V. Bernard, N. Kaiser, J. Kambor, and Ulf-G. MeiBner, Nucl. Phys. B388, 315 (1992). 

[8] S. Weinberg, Nucl. Phys. B363, 3 (1991). 

[9] S. Weinberg, Phys. Lett. B251, 288 (1990). 

[10] C. Ordonez, L. Ray and U. van Kolck, Phys. Rev. Lett. 72, 1982 (1994); Phys. Rev. 
C53, 2086 (1996). 

[11] U. van Kolck, Phys. Rev. C49, 2932 (1994). 

[12] T.-S. Park, D.-P. Min and M. Rho, Phys. Rep. 233, 341 (1993). 

[13] S. Weinberg, Phys. Lett. B295, 114 (1992). 

[14] S.R. Beane, C.Y. Lee and U. van Kolck, Phys. Rev. C52, 2914 (1995). 

[15] V. Bernard, N. Kaiser and Ulf-G. MeiBner, Z. Phys. C70, 483 (1996). 

[16] V. Bernard, N. Kaiser and Ulf-G. MeiBner, Phys. Lett. B378, 337 (1996). 

[17] V. Bernard, N. Kaiser and Ulf-G. MeiBner, Phys. Lett. B383, 116 (1996). 

[18] J.H. Koch and R.M. Woloshyn, Phys. Rev. C16, 1968 (1977). 

[19] P. Bosted and J. Laget, Nucl. Phys. A296, 413 (1978). 

[20] G. Faldt, Phys. Scrip. 22, 5 (1980); G. Faldt and T.E.O. Ericson, Phys. Lett. B89, 
173 (1980). 

[21] S.S. Kamalov, L. Tiator and C. Bennhold, Phys. Rev. Lett. 75, 1288 (1995). 

[22] R.B. Wiringa, V.G. Stoks and R. Schiavilla, Phys. Rev. C51, 38 (1995). 

20 



[23] R.V. Reid, Ann. Phys. (NY) 50, 411 (1968). 

[24] M.M. Nagels, T.A. Rijken and J.J. de Swart, Phys. Rev. D17, 768 (1978). 

[25] M. Lacombe et al., Phys. Rev. D12, 1495 (1975); Phys. Rev. C21, 861 (1980). 

[26] R. de Tourrell and D.W. Sprung, NucL Phys. A210, 193 (1973). 

[27] R. Machleidt, K. Hohnde, and Ch. Elster, Phys. Rep. 149, 1 (1987). 

[28] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

[29] P. Argan et al., Phys. Rev. Lett. 41, 629 (1978); Phys. Rev. C24, 300 (1981). 

[30] P. Argan et al., Phys. Lett. B206, 4 (1988). 

[31] M.L Levchuk, V.A. Petrun'kin and M. Schumacher, Z. Phys. A355, 317 (1996). 

[32] S.S. Kamalov, L. Tiator and C. Bennhold, Phys. Rev. C55, 98 (1997). 

[33] MAMI Proposal Al/1-96, Spokesperson: R. Neuhausen. 



Figures 




(b) 

Figure 1: Three-body interactions which contribute to neu- 
tral pion photoproduction at threshold to order q^ (in the 
Coulomb gauge). 
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(n) 
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(P) 

Figure 2: Graphs contributing at order q^ to neutral pion 
photoproduction. The hatched circles denote an insertion 
from O^fyf. The time-ordered graphs are distinguished by bold 
nucleon lines. 
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(a) 



(b) 



(c) 



Figure 3: Four-fermion terms contributing to (charged) pion 
photoproduction for z/ = (graphs a and b). Graph b is 
generated from graph c by minimal substitution. 
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Figure 4: A two-nucleon reducible relativistic Feynman 
graph. The cross denotes a generic interaction vertex. 
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(2) 



(3) 
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(4) 
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Figure 5: Decomposition of Fig. 4 into distinct time order- 
ings. Graphs (4), (5) and (6) are 1/m corrections and so the 
sum of (1), (2) and (3) is the Feynman graph in HFF. 
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Figure 6: Time shces for graphs (1) and (2) of Fig. 5 with 
pion and photon attached. All energies flow from bottom to 
top. 
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